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Abstract. By investigating wave properties at cloak boundaries, invisibility 
cloaks with arbitrary shape constructed by general coordinate transformations are 
confirmed to be perfectly invisible to the external incident wave. The differences 
between line transformed cloaks and point transformed cloaks are discussed. The 
fields in the cloak medium are found analytically to be related to the fields 
in the original space via coordinate transformation functions. At the exterior 
boundary of the cloak, it is shown that no reflection is excited even though the 
permittivity and permeability do not always have a perfect matched layer form. 
While at the inner boundary, no reflection is excited either, and in particular 
no field can penetrate into the cloaked region. However, for the inner boundary 
of any line transformed cloak, the permittivity and permeability in a specific 
tangential direction are always required to be infinitely large. Furthermore, the 
field discontinuity at the inner boundary always exists; the surface current is 
induced to make this discontinuity self-consistent. For a point transformed cloak, 
it does not experience such problems. The tangential fields at the inner boundary 
are all zero, implying no field discontinuity exists. 
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1. Introduction 

The recent exciting development of invisibility cloaks, has attracted intense attentions 
and discussions [1]- [16]. Theoretically, the cloaks are constructed easily based on 
a coordinate transformation method as proposed in Ref. [1]. The object inside 
the cloak is invisible to the outside observer, because the light is excluded from the 
object and the exterior field is not perturbed. The invisibility of the linearly radially 
transformed cylindrical and spherical cloaks has been confirmed by both numerical 
calculations [3,4] and analytical solutions [5,6]. Experimentally, the invisibility cloak 
with simplified material parameters has been implemented by Schurig et.al at the 
microwave frequency. Inspired by the idea of the invisibility cloak, some interesting 
applications, such as field concentration [9], field rotation [10], and electromagnetic 
wormholes [11], have been proposed. 

Up to now, most of discussions on invisibility cloaks focus on the cylindrical and 
spherical cloaks produced by a coordinate transformation only in the radial direction. 
For instance, in Ref. [1] , linearly radially transformed cylindrical and spherical cloaks 
are discussed in detail, and their invisibility is confirmed by ray tracing. In Refs. [5] 
and [6] , the invisibility performances of such cylindrical and spherical cloaks are further 
confirmed by obtaining the exact fields in the cloak medium directly from Maxwell's 
equations. In practice, it is sometimes desirable to have invisibility cloaks whose 
shapes are tailored for the objects to be cloaked. Thus, one needs to understand well 
about the properties of invisibility cloaks with arbitrary shape produced by general 
coordinate transformations. However, the investigations on an invisibility cloak with 
arbitrary shape are only seen in few papers [13,14] . The mechanism why the invisibility 
of a general cloak produced by compressing space is ensured, is still unclear. In this 
paper, we investigate the electromagnetic (EM) properties of invisibility cloaks with 
arbitrary shape constructed by general coordinate transformations, and we confirm 
their perfect invisibility. To figure out invisibility cloaks' main physical properties, we 
only focus on the ideal case without considering the practical implementation in this 
paper. 

The paper is organized as follows. In section II, Maxwell's equations in a curved 
coordinate system are derived. In section III, we show how to construct an invisibility 
cloak by compressing space in a general manner. In section IV, the wave behaviors 
and the medium properties at the exterior boundary of the cloak are investigated. In 
section V, we study the the wave behaviors and the medium properties at the inner 
boundary of the cloak. Through sections IV and V, the invisibility of cloaks with 
arbitrary shape is confirmed, and the fields in the cloak medium are derived with 
simple expressions. In section VI, the cloak parameters and the fields in the cloak are 
derived when the transformed space is described under arbitrary coordinate system. In 
section VII, two examples of invisibility cloaks, i.e., cylindrical and spherical invisibility 
cloaks are investigated. In section VII, the paper is summarized. 

2. Maxwell's equations in a curved coordinate system 

Maxwell's Equations in a Cartesian (x, y, z) space take the form as 

VxE=- — , VxH = — +j, V-D = p, VB = 0, (1) 

with 

D = e ? • E, B = nof ■ H. (2) 
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Consider the transformation from cartesian space to an arbitrary curved space 
described by coordinates 9i , 92 , 93 with 

x = h (t?i, 92,93), y = him, 92,93), z = 73(91,92,93)- (3) 

The length of a line element in the transformed space is given by all 2 = 
[dqi,dq2,dq3]Q[dqi,dq2,dq3] T , where the superscript "T" denotes the transpose of 
matrix, and Q — gg T with 



(4) 



The volume of a space element is expressed as dv — dct(g)dqidq 2 dq3, where det(g) 
represents the determinant of g. Here, it is noted that the way of describing the space 
transformation in this paper is similar as in Rcf. [13], where the time transformation 
is also taken into account. The space-time metric tensor g Qj a defined in Ref. [13] is 
ga/3 = diag[l, —Q] in the present paper, where only space transformation is considered. 
Then Maxwell's equations in the curved space take the form as [1, 13] 
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V,xE = --, V,xH=-+j, V ? -D = p, V 9 -B = 0(5) 



with 



D = e e • E, B = fi JI • H, (6) 

t = det( 5 )( 3 T )- 1 t 5 - 1 , % = det(g)(g T )- 1= pg-\ (7) 

j = det( 5 )( ff T )- 1 j, p = det(g)p, (8) 

E = . 9 E, H = . g H, (9) 

where the superscript "-1" denotes the inverse of matrix. 

The permittivity and permeability s and ~p in the Cartesian space are considered 
for a general case, i.e., they can be tensors. It is seen above that Maxwell's equations 
in the curved space have the same form as in the Cartesian space. However, the 
definitions of the permittivity, permeability, current density, and electric charge 
density are different, as shown in Eqs. (7) and (8). 



3. Construction of invisibility cloaks 

To construct a cloak, one usually starts from compressing an enclosed space with 
the exterior boundary unchanged [1]. As seen in Fig. 1, the region enclosed by 
boundary Si is compressed to the region bounded by the exterior boundary Si and 
the interior boundary S2. Such a space compression can be viewed as a certain 
coordinate transformation described by Eq. (3), which makes a connection between 
the points in the compressed space with coordinates (91,92,93) and the points with 
Cartesian coordinates (x, y, z) in the original space. The exterior boundary Si satisfies 
91 = x, 92 = y, 93 = z. Notice the interior boundary 5*2 is obtained by blowing up a 
line or a point [15]. Thus the cloaks can be divided into two classes: line transformed 
cloaks and point transformed cloaks. The compressed shaded region in Fig. 1 is the 
desired cloak. 

The permittivity and permeability tensors of the cloak in Cartesian coordinate 
system are given in Eq. (7). It seems that the cloak medium is very complex, 
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Figure 1. The cross section of the cloak (shaded region). 

whose permittivity and permeability are tensors and their values vary with the spatial 
location. However, the eigen functions of the wave equations in the cloak medium are 
quite simple, which relate with the eigen function of the uncompressed space by Eq. 
(9). 

In order to achieve the invisibility, the cloak should be able to exclude light from a 
protected object without perturbing the exterior field. Thus for the above cloak, at the 
exterior boundary Si , external incident light should excite no reflection. While at the 
interior boundary 5*2, no reflection is either excited, and light can't penetrate into the 
cloaked region. In the following sections, we will prove the invisibility of the cloak by 
investigating the wave behaviors at the cloak's exterior and inner boundaries. For the 
simplicity of our discussions and considering the practical application, the invisibility 
cloak is considered to be placed in air. Then the permittivity and permeability of the 
cloak in Eq. (7) will be simplified to 

l = p = det( 5 )( 5 T )-V\ (10) 

which is the same as proposed in Refs. [1] and [13]. The cloak is considered to be 
lossless at the working frequency. 

4. Cloak's exterior boundary 

In this section, we will prove that no reflection is excited at the exterior boundary. 
The transmitted electric field E 1 and magnetic field H l without interacting with the 
inner boundary S2 are expressed as 

= g~E> 1 , EP = gW, (11) 

where E l and H 1 represent the electric and magnetic fields of the external incident 
waves. According to Eq. (9), it is easily seen that the fields expressed in Eq. (11) 
satisfy Maxwell's equations in the cloak medium. Thus in order to prove no reflection 
excited at Si, one only needs to confirm that tangential components of E* (IT) and 
E 4 (EP) keep continuous across Si. 

Decompose E* and H* into E* = [E l n , E\ xl E\ 2 ] and W = [#*, W ti , H l t2 ], where 
the subscripts "n" represent Si's normal direction pointing outward from the cloak; 
"ii" and "t 2 " represent Si's two tangential directions, which are vertical with each 
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other. Thus Eq. (11) can also expressed as 



EL 



= [n, t u f^-^E 1 , 



HI 
HI 



1 2 



= [n, h, t 2 ]- 1 g^\ 



(12) 



where n, ti, t 2 represent the unit vectors in n, ti, and t 2 directions, respectively. 

At the exterior boundary Si, qi — x, q 2 — y, and q 3 = z. So fi(qi,q 2 , qs) — q% = 
(i = 1,2,3) characterize the exterior boundary Si. Therefore, it is obvious that the 
vectors V q fi~Ci (i — 1, 2, 3) lie in the same line as the normal direction n of Si, where 
Ci = x,C 2 = y, and C3 — z. For the special case when V g /, — C, = 0, V g /j — Ci can 
be expressed as On, i.e., the vector with the magnitude in the n direction. Therefore, 
g on Si can be expressed as 

g =[Fin + x, F 2 n + y, F 3 n + z\, (13) 

with 



\FA 




1)2 + ( ^L)2 + ( |/l ) 2 ) 

dqj dq k 



(14) 



where i,j, k = 1, 2, 3 and i ^ j ^ k; Fi = \Fi\ when the direction of V g /j — Ci is as 
the same as the n direction, and = —\Fi\ if the direction of V g /j — Cj is opposite 
to the n direction. Substituting Eq. (13) into Eq. (12) and noticing that n, ti, and 
t 2 are orthogonal with each other, it is easily obtained that at Si 

' .ti, Hi =U i -ti, (15) 



E l tl = E l 



El 2 =W- 1 2 , W t2 =H i ■ t 2 



(16) 



which indicates that the tangential components of E* (H 4 ) and E* (H l ) are continuous 
across Si. Thus, it is proved that no reflection is excited at the exterior boundary. 

Consider the permittivity and permeability at Si for the transformed cloak. It 
should be noticed that no flection excited at the exterior boundary does not imply that 
the exterior boundary is a perfectly matched layer (PML), where the permittivity and 
permeability at Si have the PML form s = ~p = diag[u, 1/u, 1/w] with the principle 
axes in n, ti and t 2 directions, respectively. We find that parameters at Si have the 
PML form only when g is a symmetry matrix. Observing Eq. (13), we have g T t\ = ti 
and g T t 2 = t 2 , indicating that ti and t 2 are the eigen vectors of g T with the same 
cigen value 1. Considering g is a symmetry matrix, we can know the other eigen 
vector of g T is n with eigen value det(g). Thus, n, ti, and t 2 are eigen vectors of 
Q^ 1 = {g T )~ 1 g~ 1 , with eigen values 1/ dct(.g) 2 , 1, and 1, respectively. Thus based on 

Eq. (10), e and Ji for a symmetric g can be expressed as 

t = I = dia 9i de ^ g y det(g), det(g)}, (17) 

where the diagonal elements correspond to the principle axes n, ti, and t 2 , respectively. 
The radially transformed cylindrical and spherical cloaks fall into this category [1]. 
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5. Cloak's inner boundary 

In this section, we will prove that at the inner boundary S 2 , no reflection is excited 
and no field can penetrate into the cloaked region. As discussed in section III, the 
inner boundary is constructed by blowing up a line or a point, as seen in Fig. 2(a) 
and (b). So in the following, two cases that: (1) line transformed cloaks, (2) point 
transformed cloaks, will be discussed separately. 

5.1. Case (1): line transformed cloaks 

Assume that x = b\(s) , y = b 2 (s) and z = b 3 (s) characterize the line, which is mapped 
to the inner boundary S 2 . We have fi(qi, 92, 93) = bi(s) an d /2(<7i, 92, 93) = b 2 (s), and 
/3(<Zi> 92, 93) = b 3 (s) at S 2 . Each point (61, 62, 63) on the line maps to a closed curve on 
S 2 . The parameter s can be expressed as a function of qi, 92, 93 with s — u(qi, 92, 93)- 
V 9 s = du/dqix + du/dq 2 y + du/dqs'z is the gradient of s, which points in the direction 
of the greatest increase rate of s. For \7 q bi, we have V q bi — dbi/dsV q s, where 
i = 1, 2, 3. Thus, Vqbi and V g s have the same direction. For the ease of our discussion, 




Figure 2. Illustration of cloak's inner boundary for (a) a line transformed cloak; 
(b) a point transformed cloak. 

we again decompose the incident fields at the inner boundary as E* = [E^, E% , E% ], 
H l = \H % n , , HI ] , where the subscripts "n" denotes SVs normal direction, which 
points outward from the cloaked region; "ii" and "£2" denote the tangential directions 
of ^2, with t\ vertical with the plane determined by two vectors in n and V g s 
directions, and ^ vertical with t\. Since s varies on the surface S2, the direction 
of V g s denoted by s q should not be parallel to SVs normal direction. Thus, the plane 
determined by the vectors in n and V q s directions always exists. The n, t\, and 
t 2 directions are unique, as illustrated in Fig. 2(a). E z and H l at S2 can also be 
expressed in Eq. (12), however, with the different definitions of n, t\ and t 2 . Since 
/i(9i,92,^) — bi(s) = (i = 1,2,3) characterize the inner boundary S2, V g /i — V q bi 
characterize the normal direction of S2. Then g at S2 can be written as 

g = [F x n + V„&i, F 2 n + V q b 2 , F 3 n + V q b 3 ], (18) 

with 

IF I = y/(dfi/d qi - dbi/d qi ) 2 + (df l /dq 2 - dbi/dq 2 ) 2 + (dfjdqs - dbi/dq 3 ) 2 , (19) 
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where F t = \Fi\ when the direction of V g /j — V ' q bi is as the same as n direction, and 
F, t = —\Fi\ when the direction of V g /j — V 9 &i is opposite to n direction. 

Notice that t\ is orthogonal with both n and s q , where s denotes the unit the 
vector in the direction of V g s. Substituting Eq. (18) into Eq. (12), it is easily derive 
that 

El, = H\ Y = 0. (20) 
However, the other components of fields are not zero. In particular, 

E i t2 = {s-t 2 )[B 1 ,B 2 , fl 3 ]E\ (21) 

H} 2 = (8-t 2 )[B 1 ,B 2 , B 3 }U\ (22) 

E i n = [F 1 +B 1 {s-n),F 2 + B 2 {s-n), F 3 + B 3 (s • n)]E\ (23) 

Hi = [Fi + ■ n), F 2 + B 2 (s- n), F 3 + B 3 (s ■ n)}W, (24) 

with 



B % = sjdbi/dql + dbi/dql + dh/dq 2 . (25) 

To further investigate how the waves interact with the inner boundary, the values 
of the permittivity and permeability at S 2 are needed. Observing g expressed in Eq. 
(18), it is easily obtained that g T t\ = 0, indicating that Qt\ = gg T t\ = 0. Thus one 
of Q's eigen vectors is t\ with the eigen value Xt 1 = 0, implying dct(g) = 0. Because 
Q is a symmetry matrix, the other two eigen vectors denoted by a and b should be 
orthogonal to each other and in n — t 2 plane. The corresponding eigen values are 
denoted by A a and At,, respectively, with 

A a A b = \n x s\ 2 \F x B\ 2 , (26) 

where F — F,x + F 2 y + F 3 z, B — B,x + B 2 y + B 3 z. 

Since A tl A a A b = dct(Q) = det(g) 2 , A tl = det{g) 2 / (\ a \ h ) . Observing Eq. (10), it 

is obtained that Qe = Qji = dct(g). Therefore, it is known that t\, a and b are the 

principle axes of the cloaked medium at 52 with e and ~p expressed as 

t = fJ = diag[X a X b /det(g), det(g)/X a , dct(g)/A 6 ], (27) 

where the diagonal elements correspond to the principle axes ti,a, and b, respectively. 
Since det(g) = 0, we have e a = fi a = e b = fib = 0, indicating the cloaked medium at 
S*2 is isotropic in n — t 2 plane. Therefore, n and t 2 can be considered as the principle 
axes with e n — [i n = et 2 = /it 2 = 0. For et 1 and \it x , it is seen that they have infinitely 
large values. Thus, the inner boundary operates similarly as a combination of the PEC 
(perfect electric conductor) and PMC (perfect magnetic conductor), which can support 
both electric and magnetic surface displacement currents in t\ direction [8, 12, 15]. In 
order to have zero reflection at S 2 , the boundary conditions at this PEC and PMC 
combined layer require that the incident electric (magnetic) fields in t, direction and 
normal electric (magnetic) displacement fields are all zero. From Eq. (20), we have 
E l ti = H l ti = 0. Since e„ = fi n = 0, it is obtained that D % n = B l n = 0. Therefore, it 
is achieved that no reflection is excited at 52, and E* and H l expressed in Eq. (11) 
are just the total fields in the cloak medium. The PEC and PMC combined layer 
guarantees that no field can penetrate into the cloaked region. It is worth noting 
that the induced displacement surface currents in t\ direction make E\ 2 and Hl 2 at 
52 down to zero. However, E\ 2 and H\ 2 are not zero at the location approaching 52 
in the cloak medium. Thus E\ 2 and H\ 2 are discontinuous across the inner boundary 
5 2 [8,12,15]. 
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5.2. Case (2): point transformed cloaks 

In this case, a point with the coordinate [ci,c 2 ,c 3 ] maps to the inner boundary. 
At S 2 , we have h{q\,q 2 ,q 3 ) = Ci, f 2 {qi, q 2 , 93) = c 2 , and f 3 (qi, q 2 , 93) = c 3 . The 
incident electric and magnetic Jields at^the inner boundary can be decomposed into 
E* = [E l n ,E l ti ,E t2 ] and H* = [H l n , H^, fl|J, where the definition of the subscript "n" 
denotes S 2 's normal direction, which direct outward from the cloaked region; "ti" and 
"t2 represent 6*2 's two tangential directions, which are vertical with each other, as 
shown in Fig. 2(b). Consider g at the inner boundary, which can be expressed as 

g = diag[F 1 n, F 2 n, F 3 n], (28) 

with 



I # I = ^/{dh/d qi y + {df l /d q2 y + m/d q3 r, (29) 

where i — 1,2,3. Then substituting Eq. (27) into Eq. (11), we derive that at S 2 

K = K = 0, (30) 

K = K = °> ( 31 ) 

K = [F u F 2 , F 3 ]E\ Hi = [F u F 2 , F 3 ]H\ (32) 

Unlike the case (1), in this case tangential fields are all zero, implying that no field 
discontinuity exists at S 2 . 

Analyzing Q similarly as in the case (1), we obtain that n is eigen vector of Q with 
the eigen value A„ = Ff + F% + -Ff . While the other two eigen vectors are t\ and t 2 
with the corresponding eigen values Ati = \t 2 = 0, indicating det(<?) = 0. Considering 
A„A t iA t2 = dct(g) 2 , we have \ tl = X t2 = det(g)/y/(F? + F'£ + F 3 2 ). Therefore, n, ti, 

and t 2 are principle axes of the cloak medium at S\ with e and Ti given as 



s = Jl = diag[det(g) 2 /(F{ + F 2 + F A ), ^ F 2 + F 2 + F 3 , ^ F 2 + F 2 + F% (33) 

where the diagonal elements are in the principle axes n, t\, and t 2 , respectively. Since 
det(gf) = 0, e n = fi n = 0. Considering that e„ = fi n = and tangential components 
of incident fields at S 2 are zero, it can be conclude that no reflection is excited at S 2 , 
and no field penetrates into the cloaked region. The fields expressed in Eq. (11) arc 
the total fields in the cloak medium. 

In the above sections, it has been proved that no reflection is excited at both the 
exterior boundary and the inner boundary of the cloak, and no field can penetrate 
into the cloaked region. Therefore, the invisibility of invisibility cloaks with arbitrary 
shape constructed by general coordinate transformations is confirmed. 

6. Transformation under arbitrary coordinate system 

The cloak parameters and the fields inside the cloak are expressed in Eqs. (10) and 
(11), respectively. These results are expressed under the Cartesian coordinate system, 
i.e., (91,92,93) representing Cartesian coordinates in the transformed space. However 
sometimes, it is much easier to discuss cloaks under other coordinate systems, such 
as the cylindrical cloak under the cylindrical coordinate system (r,(p,z). Thus, it is 
necessary to obtain the corresponding expressions for cloak parameters and the fields 
in a cloak under an arbitrary coordinate system, which has also been discussed in Ref. 
[13]. 
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Consider coordinate transformation, where (q lt 52,93) denotes the coordinates of 
an arbitrary coordinate system (u, v, w) in the transformed space. The spatial metric 
tensor of such coordinate system is Q u = g u g^, where g u can be obtained easily by 
considering the relationship between Cartesian coordinate system and this arbitrary 
coordinate system. For the cylindrical coordinate system and spherical coordinate 
system, the metric tensors are diag[l,r 2 , 1] and diag[l, r 2 , r 2 sin 2 9], respectively. The 
spatial metric tensor of (91,92,93) is expressed as Q q — g q g T q ', where g q is shown in 
Eq. (4). Assuming that (x\,yi,z\) are the corresponding Cartesian coordinates of 
(91,92,93) in the transformed space, the spatial metric tensor of {x\,y\,z{) is then 
obtained as Q c — g c gt, where g c = 9^\9 q , and g u \ represents g u expressed under the 
coordinates (91,92,93)- Thus, from Eq. (10), the permittivity and permeability of 
the cloak in Cartesian coordinates are obtained as e — ~p — ^(g 9 ^ 9uiQ q 1 9m- Then 
expressing e and ~p in (u, v, w) coordinate system, we easily have 



= = det(g q 
e = ijl — 



det(5„i) 
diag[pi 1 ,p 2 1 ,P3 1 } with p^ 



PiQ^QuiPr 1 , 



'Si ii 



?ul »2 



9ul t3 



(34) 

and Qui = g u i9ui- 



where Pi 

As an example, P\ — diag[l,r, 1] for the cylindrical coordinate system. 

Consider the fields in the cloak. It is easy to know that the fields expressed in 
Cartesian coordinate system are E = g c W and H = g c H l . Thus, the fields expressed 
in the (u,v,w) coordinate system are expressed as following 



E = P^gqiguofPo 1 ^ , H = PiQ- 1 1 ff,(fl u o) r i , o" 1 H < , 



(35) 



where g u o represents g u expressed under the coordinates (q 1 , q 2 , 93) and Q u q = .9«o5uo, 



P = diag[p x a , p 2 ° , _p 3 °] with p^ = 



9u0n 



9l , 



9l , 



E 1 and IT represent 



incident electrical and magnetic field vectors expressed under the (u, v,w) coordinate 
system. 

If (91,92,93) denotes the corresponding coordinates of the coordinate system 
(u, v, w) in the original space, then g q can be written as g q = g s g u o, where 



9s = 



9<J2 
dq 3 



dq 2 
dqi 

9q 2 

dq 2 

9q 2 

dq 3 



dq a 

dqi 

9q a 

dq 2 

9q 3 

dq 3 



Then Eqs. (34) and (35) can be expressed as 

det(g s ) det(g u0 ) 



e = /i = 



det(g ul ) 



Piigir'Q-^gl'QuiPr 



(36) 



(37) 



E = P^zlgsiQuvfP^ 1 , H = P 1 Q- 1 1 . 9s (Q u0 ) T Po~ 1 H J . (38) 



7. Examples: cylindrical and spherical cloaks 



In this section, based on the results obtained above, the well known radially 
transformed cylindrical and spherical cloaks will be discussed as examples. 



Coordinate transformation makes perfect invisibility cloak with arbitrary shape 10 



7.1. Cylindrical cloaks 

A two-dimensional cylindrical cloak is constructed by compressing EM fields in a 
cylindrical region r < b into a concentric cylindrical shell a < r < b. Its inner 
boundary is blown up by a straight line. Thus, a cylindrical cloak is actually a 
line transformed cloak. Here consider a generalized coordinate transformation that 
r = f(r) with f(a) = and f(b) = b, while 9 and z are kept unchanged. Thus, Q u q 
and Qui defined in the above section are diag[l, /(r) 2 , 1] and diag[l, r 2 , 1], respectively, 
which indicates that det(g„o) = f( r ) and det(g„i) = r. Po and Pi are diag[l, f(r), 1] 
and diag[l, r, 1], respectively. g s is equal to diag[f (r), 1, 1]. Substituting these 
expressions into Eq. (37), the permittivity and permeability of the cloak expressed in 
cylindrical coordinate system are obtained easily 

fir) r/'(r) f(r)f'(r) 

e r =fi r = ,,. ,e e = HB= ,,. ,e z = fi z = , (39) 

rf (r) f(r) r 

It is seen that at the exterior boundary r = b, the cloak medium has the PML form 

with e r = fi r = 1/f (b) and eg = fig = e z = fi z = f (6), which results from the 

symmetry of g c , which can be calculated easily. 

Consider the fields E and H incident upon the cloak. It is derived that 

PiQui 9s(Quo) T Po = diag[f (r), f(r)/r, 1]. Then, the fields in the cloaked medium 

can be obtained directly from Eq. (38) 

E r (r,0,z) = f'(r)Ei(f(r),0,z), H r (r,6,z) = / '(r)flj(/(r),M), ( 4 °) 
Eg{r, 9, z) = l^lEi(f(r), 6, z), Hg(r, 6, z) = Mflj(/( r ), 0, z), (41) 

E z (r,e,z) = Ei(f(r),0,z), H z (r,6,z) = flj(/(r),<M), (42) 

where \E l rl Eg,E z \ and [H^.,Hg,H z ] are components of the incident fields expressed 
in cylindrical coordinate system. When f(r) = b(r — a)/(b — a), substituting the 
expression of f(r) into Eqs. (40)-(42), we obtain the fields in the cloak medium, 
which is just the result in Rcf. [6]. 

At the inner boundary, one can easily see that s q and t^ are both in the z direction, 
ti is in the 9 direction, and n is in the r direction. Therefore, no matter what f(r) 
is, eg and fig are infinitely large, and other components are zero, as seen in Eq. (39). 
Eg and Hg, D r and B r are all zero at S*2, which guarantees that no reflection is 
excited at S2 as analyzed above. The surface displacement currents are induced to 
make E z and H z down to zero at S2. While E z and H z are not zero at the locations 
approaching S2 in the cloak medium. Thus, E z and H z are discontinuous across the 
inner boundary [8, 12, 15]. 

7.2. Spherical Cloaks 

A three-dimensional spherical cloak can be constructed by compressing EM fields in 
a spherical region r < b into a spherical shell a < r < b. Its inner boundary is 
blown up by a point. Thus, a spherical cloak is actually a point transformed cloak. 
Here a generalized radial coordinate transformation that r = f(r) with f(a) = and 
f(b) = b, is considered. Similarly as in the above example of the cylindrical cloak, the 
permittivity and permeability of the spherical cloak expressed in spherical coordinate 
system are derived 

f(r) 2 

e r = fi r = • eg = fig = e<t, = fi^ = f (r). (43) 

r z f (r) 
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At the exterior boundary r = b, the cloak medium has the PML form with e r = fi r = 
l/f'{b) and e e = \i e = e = ^ = /'(&), due to the symmetry of g c . 

Consider the incident fields E* = [E z r , E^ £j] and H' = [H l r , HI H^] incident 
upon the cloak, from Eq. (39), the fields in the cloaked medium are obtained directly 

E r (r, 6, <f) = /' (r)Ei(f(r), 6, 0), H r (r, 9, <f) = / (r)W r (f(r), 6, 0), (44) 
E e (r,6,<t>) = S -^El{f{r),e,^), H e (r,d,0) = Mflj(/( r ), 6, <f>), (45) 

Mr,9,4>) = l^Ei(f(r),0,<P), Ht(r,0,4>) = ^flj(/(r),M).(46) 

When f(r) = b(r — a)/(b — a), the fields obtained from the above equations agrees 
with the results in Ref. [5]. However, the process of the calculation here is simpler. 

At the inner boundary, observing Eqs. (45) and (46), the tangential components 
of fields axe zero [5,15,16]. Combining with. e n = n n = 0, it is known that no field can 
penetrate into the cloaked region. 

8. Conclusions 

In this paper, we have studied the properties of invisibility cloaks constructed by 
general coordinate transformations. The invisibility of cloaks is confirmed, by proving 
that no reflection is excited at both the exterior and interior boundaries of the cloak, 
and no field can penetrate into the cloaked region. The fields in the cloak medium 
are related to the fields in original EM space through E = gW and H = gH\ 
Therefore, to calculate fields in the cloak medium, there is no need to process tedious 
calculations from the complex material parameters. At the exterior boundary, when g 
is a symmetry matrix, the permittivity and permeability of the cloak medium have the 
PML form, which is just the case for our well known radially transformed cylindrical 
and spherical cloaks. At the interior boundary, the properties of the cloak for line 
and point transformed invisibility cloaks are quite different. For a line transformed 
cloak, the components of the permittivity and permeability in t\ direction (defined in 
section V) are infinitely large. While the other components are all zero. The fields 
in direction are discontinuous across the inner boundary. The surface displacement 
currents in t\ direction are induced to make this discontinuity self-consistent. For any 
point transformed cloak, at the inner boundary, the components of the permittivity 
and permeability don't have infinitely large component, and the permittivity and 
permeability in the normal direction are zero. The tangential fields at the inner 
boundary are zero, implying no discontinuity exist. Therefore, comparing to line 
transformed cloaks, point transformed cloaks are more practical due to the absence of 
the singularity of the cloak medium. 
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